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Abstract. — Non-linear sigma models that arise from the supersymmetric approach to 
disordered electron systems contain a non-compact bosonic sector. We study the model 
with target space H 2 , the two-hyperboloid with isometry group SU(1, 1), and prove that 
. in three dimensions moments of the fields are finite in the thermodynamic limit. Thus 

the non-compact symmetry SU(1, 1) is spontaneously broken. The bound on moments is 
. compatible with the presence of extended states. 
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1. Introduction 



Random-matrix ensembles such as the Gaussian Unitary Ensemble (GUE) and its 
cousins have attracted much attention in both the physics and mathematics community 
because of its many connections to statistical many-body theory, integrable systems, 
number theory and probability. This article is motivated by the study of Gaussian 
matrix ensembles which incorporate spatial structure and thus are no longer mean field 
in character. These ensembles are sometimes called band GUE models. They have 
the advantage of being mathematically more tractable than say random Schrodinger 
operators and yet they are expected to share the same qualitative features. 

About twenty-five years ago, Wegner llWllSll introduced hyperbolic non-linear sigma 
models to study band GUE models and disordered electron systems with N orbitals 
per site. In the simplest case of these sigma models the hyperbolic 'spins' indexed by 
lattice sites of Z d take values in the hyperbolic plane H 2 equipped with its SU(1, 1)- 
invariant geometry. Soon thereafter, Efetov [E| extended Wegner' s work and intro- 
duced a class of supersymmetric non-linear sigma models. 

The supersymmetric formalism has the advantage of making it possible to perform 
the disorder average and rigorously convert random matrices to a problem in statistical 
mechanics. In particular, averages of products of Green's functions become statistical 
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mechanical correlation functions. The resulting problem is how to analyse such statisti- 
cal mechanics systems. One of the main difficulties in this analysis is the non-compact 
hyperbolic symmetry identified by Wegner which naturally arises when studying spec- 
tral and transport properties of disordered systems. This paper is devoted to showing 
that a certain class of SU(1, 1) sigma models can be effectively analysed in three di- 
mensions by using horospherical coordinates and Brascamp-Lieb inequalities. 

Let A be a periodic box in Z rf (centered at 0) with volume |A| and define R(i,j) 
with i, j E A to be the elements of a Hermitian matrix drawn from the GUE. Thus the 
probability density is taken to be proportional to exp(— Tr R 2 )Y[dR(i, j). Now let J 
be a symmetric matrix with positive entries J(i,j) which are small when the distance 

— is large. Then define a band matrix H with matrix elements 

H(i,j) = y/j(i,j)R(ij). (1.1) 

If we set 

J(i,j)=(-W 2 A+iy\i,j) (1.2) 

(with A the Laplacian of the lattice A), then this corresponds to a band of width W. 
Note that J(i : j) has exponential decay e~\ l ~j\l w . Another convenient choice of J is 
given as follows. Suppose that A is tiled by identical cubes of width W. Then define 

{Jo if i and j belong to the same cube, 
J\ if i and j belong to adjacent cubes, (1.3) 
otherwise. 

In both cases the matrix elements far from the diagonal are suppressed or vanish, and 
now the dimension of the lattice plays an important role. One expects that H has qual- 
itatively the features of a random Schrodinger operator as we let the box A approach 
Z d . In fact, for the infinite- volume limit and fixed W we know that H has pure point 
spectrum (and thus localization) at all energies in one dimension, and in any dimension 
there is localization for energies in the Lifshitz tails. Large W is expected to be roughly 
inversely proportional to the strength X of the random potential. For example in one 
dimension the localization length is proportional to X for the random Schrodinger 
operator, and proportional to W 2 for a band random matrix [FM]. 

For dimension d = 3, with J given by d 1 .2ft and large W, the average local density 
of states 

p(E)=iC 1 3m((H-E-i£)~ 1 (x,x)} (1.4) 

was studied using the supersymmetric formalism in the limit when £ goes to zero. p(E) 
was shown to be smooth for E in an interval around zero and field correlations were 
proven to decay exponentially fast [D|. 

To get information about time evolution or the behavior of the eigenstates, expecta- 
tions like 01.41 ) are not sufficient. Instead one must consider 

(|(#-£-ie)- 1 (*,)>)| 2 ) (1.5) 
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and study its behavior as £ goes to zero. If (11.51) remains bounded for x = y and 
E|A|p(J?) = 1, then the eigenstates near E are extended in the sense that the L 4 norm 
of an L 2 normalized eigenstate goes to zero in the limit |A| — > oo. 

What makes (11.51 ) more difficult to analyse than (11.41) is that the absolute value 
squared eliminates oscillations and small denominators are felt when £ is small. In 
fact, even in finite volume (11.41) does not diverge as £ goes to zero. However, (11.51) di- 
verges roughly like p / where \ is the length over which the eigenfunctions are ex- 
tended. Roughly speaking, the hyperbolic symmetry emerges because the two Green's 
functions in (11.51) have energies with imaginary parts of opposite signs. Another fea- 
ture of (11.51) is that the corresponding statistical mechanics model is not expected to 
decay rapidly in 3 dimensions but rather to exhibit a Goldstone mode so that (11.51) 
should behave in the limit £ — > like l/|x — y\, the Green's function of the Laplacian 
corresponding to diffusive time evolution. 

The main purpose of this article is to analyse Wegner's non-linear sigma model (for 
one replica) as a component of Efetov's supersymmetric model. More precisely we 
study a sigma model that arises in Fyodorov's work [F|. This model is formulated on 
a lattice, whereas Wegner's model emerges upon taking a continuum limit. 

Roughly speaking supersymmetric models of disordered quantum systems have three 
sectors: the Boson-Boson, Fermion-Fermion and Boson-Fermion sectors. The B- 
B sector has the hyperbolic symmetry; this is the sector which we study. Although 
the field in this sector may potentially diverge we show that for dimension d > 3, if 
£|A|p = 1, all moments of the field remain uniformly bounded. This is the analogue (in 
the sigma model approximation) of the conjectured bound on (11.51 ) in d > 3. The F-F 
sector may also be studied in the sigma model approximation, and it corresponds to a 
classical Heisenberg model taking values in the two-dimensional sphere. Considered 
on its own, this sector has no divergence because of the compactness of the target. In 
three dimensions, the nearest neighbor Heisenberg model has an ordered state which 
may be established by using infrared bounds. The main open problem which we do 
not address in this article is the B-F sector, which couples the other two sectors in a 
supersymmetric fashion. This must be understood to obtain a complete picture of the 
SUSY models. Nevertheless, we will see that many phenomena of interest are already 
reflected in the sigma model analysed in this article. 

We now describe the hyperbolic sigma model which we shall analyse. In a periodic 
box A e 1i d (not the A of before, but the lattice of cubes that tile the original lattice), 
we associate to each lattice site jeAa matrix 

Sj W} ' • d-6) 

where Tj is subject to the conditions T*o^T = 03= diag(l, — 1) and DetT) = 1. Thus 
Tj belongs to SU(1, 1), and Sj belongs to an adjoint orbit of SU(1, 1), which may be 
identified with the symmetric space SU(1, 1)/U(1) = H 2 where U(l) is the isotropy 
subgroup generated by 1G3. 
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The action or energy of a configuration j i— > Sj is given by 

A A (S,h) = ^^(W + ^EyeA^C^) . ( L7 ) 

where £ A denotes the sum over pairs of nearest-neighbor sites of A, and h > 0. Let 
d/u(S) denote an invariant measure on SU(1, 1) /U(l) and define 

Z A (V,h) = [c-V A ^ll jeA dv(S J ). 

Expectations in this model are given by 

Z A (M) -1 / F(S)e-V A ^H jeA d M (Sj) = (F> A (M) ■ (1-8) 

Note that by M := 5o 3 = TT* = M* > 0, the SU(1, l)-orbit S = To 3 T 1 can be iden- 
tified with the positive Hermitian matrices M in SU(1, 1). Thus, TrSjSf = TrMjMj 1 
and Tra^Sj = TrA/,, and we see that A\(S,h) > 0. Since an SU(1, 1) -symmetry 
emerges at h = 0, positivity of h is needed to make the integrals exist, h corresponds 
to a magnetic field and breaks the non-compact SU(1, 1) symmetry to U(l). 
Our main result may be stated as follows. 

Theorem 1.1. — For d>3 space dimensions there is a constant Co such that 



(Tra 3 S ) 2 ) A (M)<C (1.9) 



for all (3 > 3/2 and \A\h > 1. 



Remark. — If the SU(1, 1) -symmetry were restored in the limit h — ► 0, the expecta- 
tion of the unbounded observable (Tr035o) 2 would have to diverge in that limit. Our 
result can therefore be viewed as a statement of spontaneous symmetry breaking. (A 
more detailed discussion of what it means for a non-compact symmetry to be broken 
spontaneously has recently been given in [N|.) 

Higher moments of Trc^So can also be bounded in d > 3. For d = 1, 2 we expect 
(but do not prove) that the same kind of bound holds except that we must require 

|A| " d Z\ : CIO) 

with constants C\,C2 independent of p\ 

Our proof can easily be extended to finite-range interactions for large [3; however, 
for technical reasons it does not easily extend to infinite-range interactions. 

In the next section we shall explain the relation between the band random matrices 
and the sigma models described above. Roughly speaking, the magnetic field h is 
proportional to £, the imaginary part of the energy in (11.51) . and |3 is proportional to 
W 2 , the square of the band width. In d = 1, the action (11.71 ) just describes a random 
walk on H 2 indexed by time j E 7L. 
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In the sigma-model approximation we shall see that 

H-E+ie)~ 1 (x,x) 



(H-E+ie)- l (x,x)\ 2 \DGt(H-E + ie)\ 2 



Det(#-E + ie)| 2 

(Tro 3 5 ) z ) A (p,/i). (1.11) 

Here again we identify h with £. The extra factors of the inverse determinant appear 
because we have not included the F-F and F-B sectors. 

The proof of Theorem 11.11 relies on the use of horospherical coordinates (s,t) to 
parametrize Sj E H 2 (for the details see Sect. O. The action in these coordinates is 

A A (S,h) = Y! (cosh(ti-tj) +e' i+t j(si-Sj) 2 ) +h £ (cosh^ + s 2 ^) . (1.12) 

i,jeA j€A 

The Gibbs measure now takes the form e~$ A A( s - h ) Y[jeA &tj dtjdsj. 

Note that Aa is convex in t and quadratic in s. (We mention in passing that the Bakry- 
Emery tensor Hess + Ricci for Aa is not positive in the natural hyperbolic geometry.) 
One of the key advantages of the horospherical coordinates is that we can integrate out 
the s variables thereby producing an effective action Eh(t). For (3 > 3/2 we prove that 
the Hessian of E h {t) is positive, and in fact as quadratic forms we prove 

E"(t)= Hess E h (t) > -(P- j)A+h, (1.13) 

with A the discrete Laplacian. Now the Brascamp-Lieb inequality may be applied to 
control fluctuations of the t x in terms of ( — ((3 — j)A + hj l (x,x). This is bounded in 
three dimensions provided | A\h > 1 . In one or two dimensions one must require (II . 10b . 

The remainder of this paper is organized as follows. In Sect. 12 we describe the rela- 
tion of the sigma model dl.71 ) to the random band model described by d 1 .3b following 
ideas of Fyodorov. Horospherical coordinates are introduced in Sect. |3] and convexity 
of the effective action is proved in Sect. |4j The Brascamp-Lieb inequality together 
with a Ward identity are used in Sect.|5]to obtain bounds on the t fields. The remaining 
Sects. |6] and |7] explain how to control the s field fluctuations and the /i-regularization. 

There are a number of open questions related to this paper. The primary one is to de- 
termine whether there are analogous bounds for more general hyperbolic sigma models 
such as those of higher rank. There are also problems involving averages of Green's 
functions which are not uniformly elliptic. Note that s correlations are expressed in 
terms of 

((-Ve 2t V + he')- l (x,y)) ■ (1-14) 
Since the t fields are not bounded from below, the Green's function (11.141) is not uni- 
formly elliptic. The distribution of the t fields is given by the convex effective action 
Eh. Although we obtain good bounds on (11.141) for the diagonal x = y, the off-diagonal 
bounds obtained by our methods are not sharp. 
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2. Origin of the model 

We now review how the non-linear sigma model, (11.71) and (11.81) . arises from the 
problem of computing Green's function averages for some ensembles of band random 
matrices. Aside from putting our analysis on a solid footing in random-matrix theory 
and disordered electron physics, this review will explain the origin of the hyperbolic 
target space H 2 and its Riemannian geometry. 

Readers interested only in mathematical results, not in physical motivation and back- 
ground, are invited to skip the present section; the remainder of the paper does not 
depend on it. 

2.1. Gaussian ensembles of band random matrices. — Let A c l> d be a periodic 
box as before, and assign to every site i E A one copy V/ of an ^-dimensional unitary 
vector space. (Physically speaking we are assigning N valence electron orbitals to 
every atom of a solid with hypercubic lattice structure.) The finite-dimensional Hilbert 
space V of the random-matrix model to be specified is the orthogonal sum 

V = Vi®V 2 ®...®V\ A \ . 

The basic framework we have in mind is single-electron quantum mechanics, and our 
goal is to establish control over the spectral and transport properties of certain ensem- 
bles of random Hamiltonians H. We shall take the Hamiltonians to be elements of 
Herm(V), the space of Hermitian operators on V. 

A random-matrix model is now defined by a probability distribution on Herm(V). 
Equivalently, one may specify the Fourier transform or characteristic function: 

Q(K) = (e mHK ) , (2.1) 

where (...) denotes the expectation value w.r.t. the probability distribution for the ran- 
dom Hamiltonian H. The Fourier variable is some other element K E Herm(V). 

For simplicity we shall restrict ourselves to the case of Gaussian distributions with 
zero mean, (H) = 0. If IT is the orthogonal projector on the linear subspace Vi C V, 
we take the characteristic function to be 

Q{K) = exp (4£ijLi JijTriUiKUjE^ , (2.2) 

where the coefficients Jy are real, symmetric, and non-negative (they must also be 
positive semi-definite as a quadratic form). The choice (12.21) also implies that all matrix 
entries of H are statistically independent. 

We mention in passing that the characteristic function (12.21 ) is invariant under con- 
jugation K i — > UKU 1 by U E U where U is the direct product of all the groups of 
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unitary transformations in the subspaces: 

t/=U(Vi)xU(7 2 )x...xU(VjA|). 

This means that the probability distribution for the Hamiltonian H has a local gauge 
invariance. Models of this kind were first introduced and studied by Wegner IW1I . 

By varying the lattice A, the number of orbitals N, and the variances hj, one obtains 
a large class of Hermitian random-matrix models. In particular, if d(i,j) = — 
denotes a distance function for A, and / is a rapidly decreasing positive function on 
M + of width W, the choice hj = f(\i — j\) gives an ensemble of band random matrices 
with band width W and symmetry group U = U(A^)' A L 

There exist two distinct situations where such a random-matrix ensemble is expected 
to exhibit metallic behavior (in dimension d > 3) and efficient methods of analysis are 
available. The first one occurs when the width W is large (and, say, N = I). This case 
is dealt with by the Schafer- Wegner method [S|; it will not be considered further in the 
present paper (see however [Z | for a recent review of that method). 

The second one is called the 'granular model' . There, N ^> 1 and the diagonal of 
the variance matrix hj dominates the off-diagonal: 

hi > £ hj ■ 

Each atom i E A here is to be viewed as a grain, or small metallic particle, housing 
a large number N of electron states, and the squared matrix elements for tunneling 
between grains (hj for i ^ j) are small compared to the intra-grain matrix elements 
(hi). The appropriate method to use in this case is that of Fyodorov (Sect. l23t . Metal- 
lic behavior is expected to occur when the coefficients Nhj/ y/hiJjj arQ not t0 ° small. 

Another (perhaps more familiar) way of defining the class of granular models is 
to say that one starts from matrices H drawn from the Gaussian Unitary Ensemble 
(GUE) of matrix dimension JV|A|, partitions V = C N ^ as V = V\ © V 2 © . . . © Vj A |, 
and then multiplies the variances of all matrix elements of H connecting Vj ~ C N with 
Vj ~ by the positive number hj. This is the same as the model (11.31) described in 
the introduction, with iV being equal to the volume of the cubes tiling the lattice. 

2.2. Averaging the Green's functions over the disorder. — Fixing some lattice site 
i E A, let Ag be the average absolute square of the (£,£) part of the Green's function 
with complex energy E — ie: 

A e (E,£) = (\Tr(H-E + i£)- l Tl e \ 2 } . 

Physicists have developed an approximation scheme that calculates disorder averages 
such as this one by relating them to the correlation functions of a supersymmetric non- 
linear sigma model LEI J. Here we shall address the related, but somewhat simpler 
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problem that arises from considering 

\Tr(H-E+ie)- l Il £ \ z x |Det(#-£ + ie 
|Det(#-£ + i£ 



_l„ \1 ^ r, r- . \ |-2n 

*f (E,e) = ± ^ — '- ■ (2.3) 



(n) 

Note that B) can be viewed as the average squared Green's function for a deformed 
ensemble, where the probability distribution for the Hamiltonian H is modified by 
multiplying it with In inverse powers of \Det(H — E + ie) |. Inserting these factors 
modifies the original problem, and it is far from clear how much bearing the results 
for n > 1 will have on the case n = 0. (Physicists often use a procedure called the 
replica trick, where one attempts to infer the answer for n = by analytic continuation 
from the answer for n E N. We will make no such attempt here.) However, even if the 
n > 1 situation was quite unrelated to n = 0, analysing it would still be a necessary step 
toward establishing mathematical control over the full supersymmetric theory at n = 0. 
The reason is that Efetov's supersymmetric non-linear sigma model has the effective 

(n) 

theory of B) at n = 1 for its non-compact bosonic sector. 

In order to express Brp we introduce a pair of complex fields (|) = ((]) + ,(|)_) where 
(|)± G V. The projections = and rij(|)_ = are complex Af-component 
vectors. The scalar product is given by 

If A is a matrix or linear operator acting on V with £HcA = \{A +A*) > 0, we normalize 
our Gaussian integrals over § such that 

J e -(p + ,A$+) =Der 1 (A) . 

For fixed £, £ and A define the quadratic form 

^*,A H ^) = -i(^ + ,(H-E+iE)^ + )+i^*_ y (H-E-ie)^-) . (2.4) 
Note that 9^eA > £ and the integral over (|) is therefore defined: 

Z A = J e -(PAH<?) = |Det(//-£ + i£)p 2 , 

and 

J e -W> A ^\<$> + (£)\ 2 \^(£)\ 2 = iTr^-E + i^-^^lDet^-E + iE)! -2 . (2.5) 

For general n > 1 if we set §±(j) = {^>±i(j)^±2(j), ■ ■ -A±n(j)}, this produces the 
factor \Det(H — E + ie) \~ ln and permits us to study more complicated observables 
involving several Green's functions. 

Now we can calculate the average of Za over the randomness in H by using (12.2b . 
First consider 

' ; i((^,/ty + )-i(e,ff4>_)\ = e -^LJu^( sM i sM j) ? (2.6) 
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where 

J ~\V-U)-*+U) o' (./>o (./) 

and s = 03 = diag(l, —1). For general n, 03 is replaced by the diagonal matrix s = 
diag(Id„, — Id„). Note that the In x In matrices Mj are Hermitian and non-negative; 
we say they lie in Herm + (C 2 "). 

2.3. Fyodorov's method. — Following Fyodorov [F| we choose the matrices Mj as 
our integration variables, i.e., we push forward the measure over the (]) to a measure 
over the non-negative matrices M. The push forward may be singular. However, if 
N >2n then the push forward has a density (derived in Appendix A) given by 

nl^Det^M^M,, 

where dMj denotes a (suitably normalized) Lebesgue measure on Herm+(C 2 "). Now 
set 

dfi(Mi) := Der 2 "(M,-) dMj . (2.8) 
Then we obtain Fyodorov's formula for Z\ (for general n) in the form 

Z A = J e-hWMMtiMj) J-J ^ e Tr(.,£-e)M, D ^ ^ d ^ M ^ j (2 . 9 ) 

where the integral is over the configurations {Mj} with M, > for all i e A. 

The formulation (12.91 ) offers a transparent view of the symmetries of the problem. 
Indeed, let U(n,n) be the pseudo-unitary group of complex 2n x 2n matrices T with 
inverse T 1 = sT*s. This group acts as a transformation group on the matrices Mj by 

Mj ^ TMjT* . 

Clearly, the integration domain for Mj E Herm + (C 2 ") of Fyodorov's formula ( 12.91) is 
invariant under that group action. Since |DetT| = 1 for T E U(n,n), the same holds 
true for the integration density d/u(Mj). From 

Mjs 1 — ► TMjT*s = TMjsT 1 

one sees that the function being integrated in (12.91 ) becomes invariant under the U(n, n) 
group action when the parameter £ is taken to zero. Thus U(n,n) transformations are 
global symmetries in that limit. 

In what follows we focus on the case n = 1, where the symmetry group is U(l, 1) 
or, what amounts to essentially the same, SU(1, 1). 
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2.4. The sigma-model approximation. — The exact integral representation (12.91) is 
well suited for further analysis in the granular limit which we now consider. Thus we 
now assume N ^> 1, Jij = for \i — j\ > 2, Jy = J\ > for \i — j\ = 1, and 

Ju = J > 2Jid . 

Let us first consider the diagonal terms of the action (12.91 ): 

LjeA (jJo^(sMj) 2 -iETT(sMj) -iVTrlnM,) . 
The matrices Mjs may be expressed as 

M >^'t - P l U )) T '^ ai0) 

where p\ (j), pi(j) are positive real numbers, and 7) G SU(1, 1) is determined only up 
to right multiplication by an arbitrary element in K = U(l). The measure becomes 

d^Mj) = {Mj) + 2 P2 ° f d Pl (j)dp 2 ( j) dv K (Tj) , (2.11) 

PlU) P2(J) 

where d/Jx(Tj) is an invariant measure for SU(1, 1)/U(1). The diagonal terms of the 
action can be written in terms of p\,P2'. 

\{P\ +Pi) ~iE(pi -Pi) -NQnpi +\np 2 ) . 
The critical point for E 2 < 4NJq is given by 

Pl = 27o Pn( e ) = 2Jq ' ( } 

and p2 = pi, where is the mean-field expression for the local density of states (but 
for a factor of 7t). In this approximation, 

M j s = T j ( pN ^ + const x Id. (2.13) 

1 J \ -Pn(E)J J 

The last term arising from the imaginary part of pi,p2 is of no consequence. Thus 

M jS = p N (E)Sj, (2.14) 

where Sj was defined by Sj = TjG 3 Tj~ as before. The action function of our model 
now is 

A A (S) = ±p£>(S 7 Sy) + J*£ JeA Tr (a 3 Sj) , 

with (3 = 2Ji pjf(E) and h = 2ep N (E). Similarly, by using (l2~71) and d2~T4l) . the observ- 
able appearing in (12.51) is proportional to (TrG3So) 2 m the sigma-model approximation. 
The SU(1, l)-invariant measure d/ug^Tj) is renamed to djj(Sj). 

In order to eliminate the sigma-model approximation we must control the massive 
fluctuations of pi(j),P2(j) about the saddle. Although the Gibbs measure is complex, 
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if we integrate over these eigenvalues a new effective action is produced which is real. 
This new effective action may share the desired convexity properties with Aa(S). 

3. The model in horospherical coordinates 

Having clarified the origin of the hyperbolic non-linear sigma model in disordered 
electron physics, we now begin our study of it. In the present section we introduce a 
coordinate system that takes advantage of the hyperbolic structure of H 2 and is well 
suited for the purpose of doing analysis on the sigma model, Eqs. (11.71) and (II . 8b . 

For any connected and simply connected noncompact Lie group G with semisimple 
Lie algebra there exists an Iwasawa decomposition [H] 

G = NAK , 

where K, A, and N are maximal compact, maximal Abelian and nilpotent subgroups, 
respectively. In the case at hand, namely G = SU(1, 1) with Lie algebra 

SU(1, 1) = {xiCi +X202 + 'vC303\(xi,X2,X3) G M 3 } , 

K is the U(l) subgroup generated by 1G3, and o, are the Pauli matrices. We choose 
A ~ M + to be the Abelian group generated by Oi; the nilpotent group N then is the 
one-parameter group with nilpotent generator 02 — 103. Passing to equivalence classes 
or cosets by the right action of K = U(l) on both sides of the Iwasawa decomposition, 
one gets an identification 

H 2 ~SU(l,l)/U(l)~AA-o. 

Thus the two-hyperboloid H 2 is viewed as the orbit of the one-parameter groups N and 
A acting on the coset o = K. 

Introducing two real variables s and t, we parameterize the Lie groups N and A as 

N = {n s = e*(°2- io 3)/ 2 |s Gl}, A = {a t = c tGl/2 \t G R} . 

We refer to s and t as horospherical coordinates. Their relation to the matrix S para- 
meterizing H 2 is given by 

2 2 \ 

cosh? + ye f sinh? — (is+ y)e f | 

sinhr+ (is- y)e f cosh? + ye f J 

and the SU(1, 1) -invariant metric tensor g in these coordinates takes the form 

g = jTr dS 2 = dt 2 + e 2t ds 2 . (3.2) 

How does the action of the subgroups N, A, and K on H 2 look in horospherical co- 
ordinates? (These group actions are important because they furnish global symmetries 
of the non-linear sigma model in the limit of vanishing regularization, h 0.) First of 
all, since N is a one-parameter group one has 

n So (n s a t ) ■ o = n s+So a, ■ o , 



S03 = n s a t (n s a t ) 
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son SQ G Af acts on n s a t ■ o G NA ■ o by simply translating (s,t) i— ► (s + so,t). Second, from 
the fact that 02 — 103 is an eigenvector of the commutator action [01 , •] with eigenvalue 
—2, one easily verifies 

a tQ (n s a t ) ■ o = (a tQ n s a- k) ) (a k) a t ) ■ o = n s - lQs a t+tQ ■ o , 

so a tQ G A acts by (s, t) 1— > (e _f °5, t + to). Third, the group action of K in horospherical 
coordinates is somewhat complicated and will not be considered here. 

The energy or action function of the non-linear sigma model (11.81 ) in horospherical 
coordinates is expressed by 

A a = P£ A (coshfe - tf) + - stf^i) + ^£ . £A (cosh? 7 + \s)&) . (3.3) 

The Gibbs measure is 

As expected, dfi A ,A becomes invariant under global N and A transformations, 

n s : (si , ti) f-> (si + s , ti) , 

a t : (si,ti) ^ (e" f j,-,f,- + , (3.4) 

in the limit ^ — ► 0. 

Our observable given in Theorem [TTT] may be expressed as 

(Tra 3 5 ) 2 = (2 cosh t + sle to ) 2 . (3.5) 

Note that since the action is quadratic in s, the integral over the variable sq is Gaussian 
and can be done explicitly. 



4. Integration of the s fields 

In this section we shall analyse the action (13 .31 ). Since it is quadratic in the s fields, 
they can be integrated out. Consider the interaction between the s and t fields in (13.31) 
and define 

B(s,t) =Y! A ^(si-sj) 2 = (s,Ds) A , (4.1) 

where (/,#)a = Y,iehfi8i-> an d D is a matrix corresponding to an elliptic operator with 
periodic boundary conditions and with coefficients that depend on As a quadratic 
form D is non-negative, and its matrix elements are given by 

D U = - £ Dij . (4.2) 

r-i+i 

When the variables ti all vanish, D = —A A where Aa is the discrete Laplacian of the 
lattice A with periodic boundary conditions. Although D is elliptic, it is not uniformly 
elliptic as the |f,-| may be very large. 
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Using (14.11) and integrating over the s fields we obtain an explicit expression for the 
effective action: 



Eh = P£^cosh (ti - tj?) +C h (t) + Y tj€A (-tj + hcoskt j ) > ( 4 - 3 ) 



where 



C h {t) = JlnDet (D(t) +he') +const = ~^J e'^'^^TheA^i ■ ( 4 - 4 ) 

We are going to regard the torus variables as Cartesian coordinates of IRl A l equipped 
with the canonical Euclidean geometry, and have therefore relocated the variable vol- 
ume factors e ?i from Y[e''dtj to £/,. Notice that the effective Gibbs measure 

^ E "UieAdti 

for h = is invariant under shifts ti — ► + y. This invariance is a remnant of the 
global symmetry (13.41) of the original theory, and will play an important role in later 
discussions. 

We shall first analyse a slightly different expression 

e"^ = / e-S'MSCEteA*,) n j£ A^; , (4.5) 

where the 8-function eliminates the zero mode of 5 and makes the integral exist. By 
Gaussian integration we have 

C(t) = jlnDetD(0 , 

where D > is D acting on the orthogonal complement of the constant functions. Both 
C and D depend on A and the t field, but we shall frequently omit these dependences 
for notational brevity. The effective action in t is 

Ea = $Y! A cosh(ti-t f ) +C(t)+Y, jeA (-tj + hcoshtj) . (4.6) 

We shall set the factor of (3/2 appearing in (14.51) equal to 1. By scaling in s this simply 
shifts C(t) by a trivial constant. 

For a function F of ti (i E A) let the Euclidean Hessian of F be denoted by F": 

d 2 F 
dtjdtj 

Theorem 4.1. — For any value of the coupling parameter (3 > 3/2 and dimension 
d>\ the function E\ is convex and 

E'jl>-($-±)A A + h>-A A + h. 
Proof. — Clearly from (14.61) 

E' A > -$A A + h + C" , (4.7) 
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so it suffices to estimate C" . From (14.51) with (3/2 set to 1 we have 



(4.8) 



where (■) denotes the average over the s field with Gibbs weight e s 8. For — > 1 



we have 



and for \i — j\ = 1 



a 2 c 

dtidtj 
dtidtj 



{U i \U J )={U i U J )-{Ut)tU J ) l 



(Ur,Uj)-(^( Si - SJ r) 



while on the diagonal 



Let 



32C 



(Uf,Ut)-(U t ),. 



(4.9) 
(4.10) 



K ij = {U i ;U j ) s , 
and decompose C" into two pieces: 

^=(2<tf,) f 8y-Zy)+*y, 

where R is a local remainder term. 

By explicit computation Kij(t) > for all i, j and field configurations t, because the 
square of a Green's function arises. 

Lemma 4.2. — With U and K defined by (14. #D an J ( 14.91 ) we have for every t 

LjKii 2(U i)s . (4.11) 
TTu's relation implies that as a quadratic form 

2(t/,) s 8, ; -£, ; >0. (4.12) 

Proof. — To get (14.1 II) make the change of variables st — > (k E A). Then {/,■ — > 
f/ ( e 2Y and 5 — > 5e y , while the expectation value ([/,•) remains invariant. Differentiat- 
ing (£/,•) with respect to y at y = yields 

= 2(Ui) g -2(Ui;B), 

and since 5 = ^ £ / ^0 we obtain ( 14.1 II ). 

The non-negativity of the quadratic form (14.121 ) now follows from the Schwarz in- 
equality: 

EWi < \L K ijff\ 4 |l%/f * = 2l^./f • 
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Here we used the pointwise positivity of Kjj to write the first expression as a scalar 
product of two vectors u and v with components Ujj = y/Ki]fj and Vy = ^Kijf. □ 

We now must estimate the remaining local part R = C" — 2(U)8 + K, which is ex- 
pressed by 

R tj = -MSij + ^Jdsi-SjftSdi-jl - 1) . 
Note that R < as a quadratic form and that J^jRij = for each i. 

Lemma 4.3. — For all real f 

Y,*i]fif]\<hYFf)h (4-13) 
where V/ denotes the discrete gradient of the lattice A. 

Proof. — The left-hand side of (14.131) can be written as a sum over nearest-neighbor 
pairs ij: 

H\i-j\ = l Ri j(fi-fj) 2 ■ 

It therefore suffices to show that for each pair j we have 

^(( Si - Sj ) 2 ) s <l/2. 

This result follows from the fact that ( • ) s is a Gaussian expectation in s with terms 
e' i+t J (si — Sj) 2 appearing in the action, B. Indeed, if w, v\,V2,- ■■ are real variables it is 
a general fact that 

fcU 2 0-"' 2 -Q^ V UuUadVa ^ 
f e -cu 2 -Q(u,v) du Yl a dv a 

for any positive constant c and any Q > which is quadratic in w,v. If we set cu 2 = 
(si — Sj) e ti+tj and Q = B — cu 2 (restricted to the linear subspace given by the constraint 
E/eA^i = 0)' we obtain the desired result. □ 

From the decomposition (14.101) and the two lemmas we have 

C">R> ^A A . 

Inserting this bound into (14.71) completes the proof of Theorem 14. II 



5. The Brascamp-Lieb Inequality 

We now state the Brascamp-Lieb inequality [BJ in a form in which we shall apply 
it. Let A = A[t) be a convex function of N variables t = (ti,...,tff) E M. N , where M. N is 
the Euclidean vector space with scalar product 



(<p»o=E f= i<p<^ 
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With the function A associate the measure d/JA(t) = e A ^> Yldti- Assume that the Eu- 
clidean Hessian of A satisfies 

A"(t)>H>0, (5.1) 
where H is a positive N x N matrix independent of t. 

Theorem 5.1 (Brascamp-Lieb). — If A satisfies rt5.il) then 

{w} _ f&'WWdtj {{W])A i ({p ,„- 1(p) 

For our application we identify A with the function E = E\(t) given by (14.61) . The- 
orem |47T] tells us to put H = — (|3 — ^)A A + If we then set 

Gy=(-(P-±)A A +fc) -1 .. (5.3) 

and fix a site i E A with field variable ti, we have 

( e ctt,-a(f;> £ ^ < e ^a 2 G !7 _ (5 4) 

In dimension J > 3 and for [3 > 3/2 the Green's function Gu is uniformly bounded as 
A provided h | A| > 1. 

Now we drop the subscript E and let (■) = (-)e- Large fluctuations of the field t 
away from its average are very unlikely: 

_ p 2 

p = Prob £ {^ - (ti) > p} < e *5s . 
Indeed, when ^ — (t{) > p we have e a ( fi_ ^ _p ) > 1 for all a > 0, so that 

p < ^-('i)-p)) < e -ap e i« 2 G i; = e -4j 

for a = p/Gu. The same estimate applies to the probability of an event ti — (ti) < — p, 
so altogether we have 

Prob £ {k,-(^>|>p}<2e ^ . (5.5) 

Our estimates on the t field will be complete once we have estimated the average 
(ti). To do this consider the change of variables tj — > fy+y, Sj — » e^Sj (j G A). Then 
if we take the derivative in y of the logarithm of the partition function at y = we get 

fc£ jeA (sinh* 7 ) = 1 , (5.6) 

where 1 is produced from the 8-function: 8(e~ Y £.Sj) = e^dC^Sj). By translation in- 
variance we see that for h\A\ = 1 we have (sinh£ ( ) = 1, so that by Jensen's inequality 
and (15.41) we have 

e (f,-) < ^ = 2 + ( e - ?i ) < 2 + e"^e5 G " . 
This gives an upper bound to (tj) : 

(ti) < 1 + \G U ■ 
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To obtain the lower bound we use (sinhf,-) > 0: 

e - fo> < (e~' f ) < (e f '') <e^e2 G " . 

Hence 

-\Gu<{ti) <\ + \G H . (5.7) 
This completes our estimates on the t fields and its fluctuations. 



6. Bounds on the s fields. 

Recall that in addition to the t variables the observable given by (13.51) contains factors 
of Sq and Sq. These averages may be explicitly calculated in terms of the co variance 
for the s field given, see ( 14. II) . by D 1-1 on the orthogonal complement of the constant 
functions which we have denoted by D~ l . Here we show how to deal with the Sq term: 

(4) J = 35- 1 (0,0) 2 . (6.1) 

The Sq term is similar and can be handled in the same way. 

The operator D is non-negative but depends on t. If all tj > then D > — A a and we 
have bounds on D 1 in terms of the free Green's function — A A l in dimension three. 
However, the t field may take large negative values and so there is no uniform bound 
on D~ l (t). The control of (D _1 (0,0) 2 ) will come from the fact that large negative 
values of t are very rare by (15.51) and (I5.7I) . 

To bound the average of D _1 (0,0) 2 we shall first consider an elliptic operator L 
whose quadratic form is 

where a j > (1 + Let L denote the projection on the orthogonal complement of 

the constant functions. 

Lemma 6.1. — For d>3 and p < d — 2 the Green 's function ofL is uniformly bounded 
(as A ^Z d ) by O^L- 1 (0,0) <A p <oo. 

Proof. — Let C„ denote the cube of side 2 n centered at the origin and let % n be its 
indicator function. Note that 

fn = 2- dn Xn-2- d ^+ l h n+l 

has zero average and the square of its L 2 norm is bounded by 2~ dn . 

As a quadratic form, L restricted to C n (with Neumann boundary conditions) is at 
least 2~( 2+p )", and we therefore have 

(/„_! ,Z-Vn-l) < 2^P)-2~ d ^ . 



18 



T. SPENCER & M.R. ZIRNBAUER 



To complete the proof of the lemma note that the projection of 80 onto the orthogonal 
complement of the constants can be written as a sum over the f n . By the Schwarz 
inequality we have 

E -1 (0,0) < (ir=o(/«' Z_1 /n) 1/2 ) 2 <A P , 
provided that p < d — 2. □ 



Lemma 6.2. — There is a constant cq so that 

(p-\W) 2 )<co. (6.2) 

Proof. — Fix some value of p with < p < j, and for each integer k let %k(t) denote 
the characteristic function of the set of configurations t = {tj}jeA that satisfy 

e'* + ^>e-*(|jc| + l)- p 

for all nearest neighbors x,y. We then have Xk(t)D~ l (0,0) 2 < (Ae ) 2 by Lemma |6~T1 
We now claim that for all k > K = max{0, — (t x + t y )}, 

(l-% k (t)) t <Bo- c ^ 2 , (6.3) 
which by Borel-Cantelli gives the desired statement: 

(D-^O) 2 } = ^( %0 + %1 (l- Xo ) +X2 (l- %1 ) + ...)^ 1 (0,0) 2 ) 
< const +A 2 B ^ e 2(k+i) e -c(k-K)i < const 

k>K 

To establish our claim (16.31 ) suppose that 

p -k 

J x +t y <■ c <■ e -{k+np) 

-(\x\ + l)P- 

for x in the range e' 1 < \x\ < e n+1 . Then by (I5.5I) and (I5.7I) the probability of this event 
for k > K is less than 

Q (n+l)d e -c(k-K+np) 2 

whose sum over n is no greater than e - c ( fc - K ) 2 times a constant B. □ 

The lemma works for any power of D -1 (0,0). Thus by the Schwarz inequality and 
(I5.4I) we can bound the expectation of our observable e 2t °SQ. 
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7. Adjusting the regularization 

We have used the 8-function regularization in the s variables rather than the correct 
term /?Lj £ A^e^ which appears in the action A^. 

Recall C/ 7 = ^lnDet(D + /ie / ) and C = ^lnDetD. We shall express Q in terms of 
C. To do this, let Pq denote the orthogonal projector on the vector space spanned by 
the normalized constant function \]/oj = |A| -1 / 2 , and let P = 1 — Pq. The determinant 
can be calculated in terms of P and Pq blocks: 

Det(D + /*e f ) = Vet(D + hP t ) Oo,W\|/ ) , 

where P t is given by 

P t = Pe l P -P^Pq^P • (\|/ ,eVo) _1 - 

Using the Schwarz inequality it is easy to see that P t > 0. 

Thus C h = C+ ±Trln(l +hD^ 1 P t ) + ±ln (\|/o,e f \|/ ). We have left out the ln/z term 
since it is cancelled in the normalization. 

Let F(t) be our (positive) observable. Now we can write 

(F(t)) _ (F(t)e-«») E 
{F{t))E h - {e _ Rh)E , 

where 

R h = C h -C= iTrln^+^D-^ + ilnCxi/o.eVo) • 
Since the first term of Rf, is positive we have 

< ^(0(¥o ; e> )^) £ < {F(t)e -^o) )E e (R h ) E 

where we have used Jensen's inequality. Since P t is positive we have 

(Tr\n(l +hD- 1 P t )) E < h(Tr(D- 1 P t )) E < h\A\(e to D-\0,0)) E < const, 

and the other term in (Rh) is estimated by 

(iln(\]A ,efyo)) < ±(e'°) < const . 

The desired bound on (F)^, now follows from estimates we obtained for ( This 
completes our proof of Theorem ll.il 



8. Appendix: Push forward of measure in Fyodorov's method 

Consider the mapping 

y : Hom(C' 7 ,C w ) -> Herm+(C' 1 ) , (p ^ <p*<p = M , 
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and fix some (translation-invariant) Lebesgue measure dq>dq> on Hom(C B ,C iV ). We 
claim that, if N > n, there exists a Lebesgue measure c„,n dM = dM n ^ (with normal- 
ization constant depending on n and N) such that the equality 

J F((p»J(pJ(p = J F(M)D&{ N - n (M)dM n ^ N (8.1) 

Hom(C",C A ') Herm+(C") 

holds for all functions M i-> F(Af) on Herm + (C") with finite integral /F(<p*<p) dcpdcp. 
In other words, \|/ pushes the measure dcpdcp forward to 

\\f(dqdip) = Det Ar_n (M) rfM M;iV . 

While this claim can be viewed and proved as a statement in invariant theory, the most 
elementary proof is to express the integrals on both sides in terms of generalized polar 
coordinates, as follows. 

Given any complex rectangular matrix cp G Hom(C",C A ') for N > n, consider the 
non-negative Hermitian matrices M = (p*(p and M' = (pep*, which are of size nxn and 
N xN respectively. The rank of M' cannot exceed n, so there must be at least N — n 
zero eigenvalues. The other n eigenvalues are in general non-zero, and coincide with 
the eigenvalues of M = cp*(p. Denote these eigenvalues by A4, . . . ,A, M ; their positive 
square roots are sometimes called the singular values of (p. There always exist 
two unitary matrices U G U(rc) and V G U(iV) such that 

9* = uV^V- 1 , (p = Vv / Xc/" 1 , 

where is the rectangular N x n matrix with diagonal entries \^ki, \f 7 k n and 
zeroes everywhere else. 

Let J(VX) be the Jacobian of this singular value (or polar) decomposition: 

j(^)= n (v^i-Vh')\v^i+Vh') 2 fiVh 1+2{N ~ n) . 

l<i<i'<n k=l 

Fix the values of n and N > n. Then by a standard argument there exists some (fixed) 
choice of Haar measure dU for U(n) such that 



/ 



F((p»J(pJ(p= / (/ F(UW- l )dU )J(VX)T\ d\f\ k 

JRI \Ju{n) J t=\ 



Hom(C",C A ') 

holds for all integrable F. Here X = diag(Xi , A,2, . . . , A, n ). 

On the other hand, fix some Lebesgue measure dM for Herm(C") . By diagonalizing 
the Hermitian matrix M and transforming the integral / f(M)dM to the eigenvalue 
representation M = UXU~ l you get 

/ f(M)dM=b n . N f (f nuw-^du^iKh-ij^Yidh. 

Herm+(C») + 3 
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The constant b n ^j is determined by the (arbitrary) choice of dM relative to dU. Now 
put f{M) = Det N n (M)F(M). Since Det^'^M) = Uk^~ n , the desired statement 
(18.11) follows (with c„_n = 2~ n /b n _N) by comparing expressions and noting 

The relation (18.11) can also be viewed from another perspective, which we shall now 
offer. First note that the integral on the left-hand side can be regarded as a distribution 
(or continuous linear functional), say /j, on Herm + (C"): 

|i:Fn J F(<p*<p) dydy. 

Next observe that the non-compact Lie group GL(n, C) acts transitively on the positive 
Hermitian n x n matrices M by 

M i— > TMT* (T G GL(n,C)) . 

Via this action we can identify (a dense open subset of) Herm + (C") with the non- 
compact symmetric space GL(n,C)/U(«). The corresponding action on functions, 
F i— > T F, is given by 

T F(M) = F(T~ l MT~ 1 *) . 
Given that T E GL(n,C) acts on the functions, there is also an action fj i— > T(jj) on 
the distributions, by T(fi) [F] = /u[ T F] . Since the Jacobian of the transformation cp i— ► 
cp o T* and cp* i-> T o cp* is Det N (T*T), the distribution /j at hand satisfies 

T(fi) =Det N (T*T)p. 

Now we make the same considerations on the right-hand side of (18.11) . i.e. for the 
distribution 

^:Fh J F(M) Det Ar_n (M) dM . 

Under the transformation M i— > TMT* the Lebesgue measure dM transforms into 
Det"(r*r) dM. Hence Det~"(M) dM is invariant under such transformations, and 

T(Ji) =Det w (rr)£, 

i.e., Ji transforms in exactly the same way as fj. It is an invariant-theoretic fact — 
resulting from the interpretation of /j as a linear functional on U(N) -invariant state 
vectors F in the oscillator or Shale-Weil representation of Sp(2nN) — that the vector 
space of distributions with this transformation property has dimension one. Therefore, 
there exists some constant c n # such that 

V = c n ,N x V ■ 

As a corollary, we note that Det~" (M) dM is an invariant measure for the symmetric 
space of positive Hermitian n x n matrices, Herm + (C") ~ GL(«,C) /U(n). The case 
encountered in the main text of the paper is obtained by replacing n — > In. 
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